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Abstract 

The data set collected with the ALEPH detector from 1991 to 1995 at LEP has been 
analysed to measure the charm forward-backward asymmetry at the Z. Out of a total 
of 4.1 million hadronic Z decays, about 36000 high momentum D* + , D + and D° decays 
were reconstructed, of which 80% originate from Z —>■ cc events. The forward-backward 
asymmetry was measured at three energy points: 

A C FB (^ = 89.37 GeV) = (-1.0 ±4.4)% 
A FB (^ = 91.22 GeV) = (6.3 ±1.0)% 
A c fb(Vs = 92.96 GeV) = (11.0 ± 3.4)% . 

From this analysis, a value of the effective electroweak mixing angle sin 2 9^ = 
0.2321 ± 0.0016 is extracted. 
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1 Introduction 



The forward-backward asymmetry in Z — > cc decays provides a direct and precise test 
of the coupling of the Z to up-type quarks. The asymmetry A FB in Z — > ff decays 
arises from parity violation in Z production and decay. In the Standard Model, the 
differential cross section, expressed as a function of the angle 8 between the outgoing 
fermion and the incoming electron, is 

--^- = 3 -(i + cos 2 e) + A { FB cose (i) 

a a cos t) 8 

At the Z pole, for unpolarised e + e~ beams, A FB is related to the pole asymmetry A° FB , 
defined in terms of the effective couplings in the improved Born approximation as 

,f 3 2g Ve gAe 2g V fg A f . s 

FB " 4 (g^ e + gL) (g vf + g| f ) • 1 J 

The measurement can be interpreted in terms of the effective electroweak mixing angle 
sin 2 ^ = i(l-g Ve /g Ae ). 

The Z — > cc decays provide a convenient way to directly test the Z coupling to up- 
type quarks, since the Z — > uu events are much harder to isolate. In this paper a sample 
of Z — > cc decays is selected using high energy D + , D° and D* + , fully reconstructed 
from their decays products. The restricted number of usable decay channels limits 
the tagging efficiency On the other hand, the high purity obtained in the selection, 
together with a minor dependence on the decay models, allow a measurement with 
small systematic uncertainties. 

First the selection of charmed mesons is described, then the background estimate 
is discussed and the asymmetry is measured on the selected sample at three different 
centre of mass energies: at the Z peak and at ±2 GeV off-peak. 



2 Reconstruction of D^*) mesons 

A detailed description of the ALEPH detector and its performance can be found in 
Ref. 0, |J. Charged particles are detected in the central part of the detector, consisting 
of a two-layer silicon vertex detector with double-sided (r-0 and z) readout, a cylindrical 
drift chamber and a large time projection chamber (tpc), which together measure up to 
33 coordinates along the charged particle trajectories. Tracking is performed in a 1.5 T 
axial magnetic field provided by a superconducting solenoid. The tpc also provides 
up to 338 measurements of ionization (dE/dx) allowing particle identification. The 
electromagnetic calorimeter is a lead/ wire- chamber sandwich operated in proportional 
mode. It is read out in projective towers of typically 15 x 15 mrad 2 size segmented in 
three longitudinal sections. The iron return yoke is instrumented with streamer tubes to 
provide a measurement of the hadronic energy. An energy flow algorithm || combines 
charged particles momenta and calorimetric energy measurements and provides a list 
of energy flow particles on which the analysis is based. 

About four million hadronic Z decays are selected as described in Ref. ||, out of the 
data set collected by ALEPH during the 1991-1995 running period at the Z resonance. 
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Table 1: Results of the D + ; D° and D* + reconstruction. The second column shows 
the number of candidates in the different channels, while the last column shows the 
number of signal events, after combinatorial background subtraction, together with the 
statistical and systematic uncertainties. 



Decay channel 


Candidates 


Signal 


D*+ -> 7T+K-7T+ 


5022 


4434± 71± 55 


D*+ -> TT+K-TT+TT 


7327 


5429± 86±124 


D° -> K-7T+ 


7682 


5032± 97± 75 


D*+ -> 7r s +K-7r+7r+7r- 


14565 


8710±121±276 


D*+ -> 7T s + K-7r + (7r°) 


10766 


5824±104±322 


D+ -> K-7T + 7T + 


12664 


6357±113±102 



Charmed mesons are reconstructed in the following decay modes (charge-conjugate 
modes are implied throughout the paper): 

(i) D*+ -> D 7T+ (ii) D*+ -> D°7T+ 

UK-7T+ UK-7T + 7T 

(ziz) D° -> K-7T+ (iu) D*+ -> D°7T S + 

'— >K~7r + 7r + 7r~ 



u) D*+ -> D°7T+ (w) D+ -> K" 

UK-7T++ (7T°) 



7T 1 71" 



where the soft pion from D* + decay is indicated as 7r+. Channel (v) is selected without 
reconstructing the 7r° by using the kinematic properties of the underlying resonances 
in the D° decay that make the K~ir + invariant mass peak near 1.6 GeV/c 2 . The 
reconstruction proceeds from channel (i) to channel (vi); once a candidate has been 
found in one event, this event is excluded in the following channels. Consequently, 
the D° — > K~7T + decays from D* + are removed from the inclusive sample (iii), thus 
avoiding double counting. The number of candidates for each selected sample is listed 
in Table 0. 

All combinations of two and four tracks or two tracks and a tt°, with total null 
charge, are considered as D° candidates, and all combinations of three tracks with 
total charge +1 are considered as D + candidates. The ir° candidates are selected from 
two-photon combinations having a x 2 probability of at least 5% for a mass-constrained 
kinematical fit 0. The invariant mass of the D candidates, with mass assignment 
according to particle hypothesis, is required to be close to the D meson mass within two 
times the invariant mass resolution. The D° — > K~tt + + (tt°) channel is reconstructed 
as the D° — > K~7r + decays, except that candidates are kept if the K~ir + invariant mass 
is between 1.5 GeV/c 2 and 1.7 GeV/c 2 . 

D* + candidates from D* + — > 7T+D decays are selected by adding an extra track 
with momentum less than 3.5 GeV/c to a D° candidate. In the channels (ii), (iv) and 
(v) the combinatorial background is reduced by requiring the D° candidate to satisfy 
| cos#*| < 0.8, where 9*, in the D° rest frame, is the angle between the D° direction 
and the sphericity axis of the three (D° — > K~7r + 7r°) or four (D° — > K _ 7r + 7r + 7r _ ) decay 
products, or the kaon direction in the decay D° — > K^7i + (n°) with an undetected n°. 
The 7r+ momentum is required to be greater than 1.5 GeV/c, favoring high momentum 
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D* +, s, in order to reject combinatorial background and Z — > bb events in which a 
b hadron decays into a D* + . In Figure [l] the mass difference AM = M D »+ — M D o 
distributions are shown. Candidates are selected in the AM region 143.5 MeV/c 2 to 
147.5 MeV/c 2 for the channels (i), (ii) and (iv) and 141 MeV/c 2 to 152 MeV/c 2 for 
channel (v). 

Candidates from D° — > K~7r + decays, channel (Hi), are selected if the kaon track 
momentum is greater than 2.5 GeV/c and the pion track momentum greater than 
1.5 GeV/c. A common vertex is searched for and candidates are kept if a vertex 
with a x 2 probability greater than 1% is found, and the projected decay length 
significance is greater than unity. The asymmetry measurement in this channel suffers 
from the presence of a significant contribution of fake candidates due to incorrect 
mass assignments, which reverse the charge assignment of the charm quark. For a 
large fraction of these candidates, the correct mass assignment is also selected in the 
event. The dE/dx measurements of the two tracks are used to choose between the two 
mass combinations. The probability to be a kaon (Pk) or a pion (P n ) is computed 
from the measured track ionization and the expectation for a kaon or a pion. The 
mass assignment which gives the highest probability Pk x P-n is kept. If no dE/dx 
measurement is available for the tracks, the choice is made randomly. This criterion 
reduces the contribution from incorrect mass assignments to 4% in the signal region. 

D + candidates, channel (vi), are selected if the kaon track has a momentum greater 
than 2.5 GeV/c and if the dE/dx measurement is more consistent with the expectation 
for a kaon than for a pion. One of the two pion tracks is required to have a momentum 
greater than 1.5 GeV/c, the other pion momentum being greater than 0.75 GeV/c. 
The three tracks are required to form a common vertex with a \ 2 probability greater 
than 1%, and a significance of the decay length, projected along the D + momentum, 
greater than 1.5. Finally, in case of multiple candidates, only the candidate with the 
largest decay length significance is kept. 

The contribution of the Z — > cc process to the D° — > K~7r + and D + — > K~7r + 7r + 
signals is enhanced to 80% by selecting D° and D + candidates with energy greater than 
half the beam energy. The resulting invariant mass distributions are shown in Figure ^| 
for the D° — > K~7r + sample and in Figure |3| for the D + — > K~7r + 7r + sample. 

3 Combinatorial background estimate 

The fraction of combinatorial background events in the D* + sample is estimated 
from the mass difference distributions. The data sample contains, in addition to 
D* +, s which are correctly reconstructed, a combinatorial background and a fraction 
of D* +, s obtained from soft pions and partially reconstructed or fake D°'s. The 
latter contribution, clearly seen in Figure 1(c), carries the correct charm quark charge, 
so it is treated as signal. The mass difference of the combinatorial background is 
obtained from D° candidates in Monte Carlo simulated events in which no D* + 's have 
been produced. A track from fragmentation is added to such candidates and the 
combinatorial background is estimated from the resulting AM distribution, normalized 
to the data in the region AM > 0.16 GeV/c 2 . In the background normalisation 
procedure reflections of the signal in the AM > 0.16 GeV/c 2 region have to be taken 
into account. This is done by using, both in data and in background Monte Carlo, 
only events in which a signal candidate is not found. 
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The fraction of combinatorial background events in the D° and D + samples is 
extracted from a fit to the invariant mass distributions (Figures [| and |3|). The 
D° and D + signals are parametrized by two Gaussians with a common mean and 
the combinatorial background by a polynomial function. Resonant background 
contributions, such as D° — > K~K + , 7r + 7r~(7r°) and D° — > K~7r + where the two mass 
assignments are reversed (D° channel), — > (frn^ and — > K*K (D + channel), are 
taken into account in the fit. Their shapes and sizes are fixed by the Monte Carlo 
simulation with branching ratio according to PDG values |4]]. 

The fitted numbers of signal events in the different channels are shown in Table [|. 

4 Measurement of the forward- backward asymmetry 

The measurement of the differential cross section for Z — > cc events (Eq. |I|) requires the 
evaluation of the angle 9 between the charm quark direction and the incident electron 
beam. This is measured from the thrust axis, oriented along the candidate direction, 
cos 9 = — Q cos^thmst, Q being the electric charge of the reconstructed K in the D° or 
D + decays. 

Together with charm events two possible sources contribute to the observed 
asymmetry: combinatorial background and charm mesons from Z — > bb events. 
Therefore the observed asymmetry in the selected sample is 

Ap£j — fsigfcAp B + /sig(l — fc)Ap B + (1 — fsig)Ap B 

where / s ; g is the fraction of D mesons in the sample, f c is the fraction originating 
from direct charm production and A h p B is the forward-backward asymmetry of the 
combinatorial background. 

The fraction, f c , of D mesons originating from direct charm production is measured 
directly from data ||. The event is divided into two hemispheres according to the 
thrust axis. A lifetime-mass tag |J is applied on the hemisphere opposite to the D 
meson, to select b hemispheres with 99% purity. The fraction, /b-tag, of D mesons that 

r 

survive the b-tag cut is used to extract the charm fraction f c = ^ j^tag w h ere 

ta ^ ebb - e c e ' 

e b b and e C c are the b-tag efficiency for b and charm events; e cE is obtained from Monte 

Carlo simulation, while e bb - is measured from an unbiased Z — > bb data sample || and 

corrected to take into account the presence of an energetic D meson in the opposite 

hemisphere. The systematic error on f c arises mainly from the uncertainty on this 

correction ||. The measured charm fractions are listed in Table |2|. 

The b asymmetry A FB is fixed to the Standard Model values at the three centre of 

mass energies, as listed in Table [| together with the dependence of the fitted charm 

asymmetry. The effective asymmetry that enters in this analysis is diluted by a factor 

(1 — 2x mix ), due to the mixing of neutral b mesons. The mixing probability Xmix 

is different for each reconstructed D meson, since it depends on the fraction of such 

mesons produced in B° decays among all b — > D* 1 * 1 , D°, D + decays. These fractions 

are derived from Monte Carlo simulations and the associated errors are taken as the 

difference between the Monte Carlo prediction and an estimate based on experimental 

measurements in the semileptonic sector 0. These values, together with the world 

average value of Xd, are used to obtain the Xmix values, shown in Table [| The 

contribution from the double-charm decays of b hadrons, in which the charge of 
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Table 2: Charm fraction and mixing probability used to extract the charm asymmetry 
from the observed asymmetry. The first column shows the charm fractions, as measured 
from data, with the statistical and systematic uncertainties. The second column shows 
the obtained fractions of D*^ 's, D° 's and D + 's coming from B° decays among all 
b — > D* ± ,D°,D + decays. The resulting values o/Xmix ore listed in the last column. 



Decay channel 


/c 


B°^D 
b^D 


Xmix 


D*+ 


— > 7T + K 7T + 


0.741 ±0.019 ±0.007 


0.80 ±0.05 


0.16 ±0.04 


D*+ 


— > TT + K~7T + 7T° 


0.743 ±0.019 ±0.007 


0.80 ±0.05 


0.16 ±0.04 


D° - 


-> K-7T + 


0.787 ±0.019 ±0.006 


0.22 ±0.05 


0.035 ±0.005 


D*+ 


— > 7T + K~7T + 7T + 7T~ 


0.783 ±0.016 ±0.006 


0.80 ±0.05 


0.16 ±0.04 


D*+ 


-> 7r+K-7r+(7r°) 


0.766 ±0.021 ±0.007 


0.80 ±0.05 


0.16 ±0.04 


D+ - 


-> K~7T + 7T + 


0.797 ±0.020 ±0.006 


0.72 ±0.03 


0.12 ±0.02 



Table 3: Value of A FB at the three energy points, used to extract A C FB . These value 
correspond to the Standard Model prediction with m t = 175 GeV/c 2 ; m# = 127 GeV/c 2 
and a s = 0.120, without QCD correction in the final state (a discussion of the 
corrections to the asymmetry in this analysis follows in the text). The last column 
shows the dependence of the charm asymmetry on the value of A FB . 



v^(GeV) 


A b FB (%) 


dAp B / dAp B 


89.37 


5.7 


-0.22 


91.22 


9.7 


-0.22 


92.96 


12.1 


-0.22 



the reconstructed has the opposite sign with respect to single-charm decays, is 
negligible due to the low momentum of the decay products ||. 

The doubly Cabibbo suppressed decay D° — > K + 7r~, which affects both the charm 
and the b component of the inclusive D° sample, has a negligible effect on the results. 
The asymmetries are corrected to take into account the fraction of the selected D'*''s 
that originate from gluon splitting, estimated to be (0.8 ± 0.4)% M. 

The asymmetry of the combinatorial background is measured from the upper side 
band of the mass peaks. Within the side bands, multiply-counted events are selected 
only once by choosing the candidate randomly. In the case of the D + , the contributions 
of the resonant backgrounds are negligible in the upper side band region. On the 
other hand, in the case of the D° — > K~tt + channel, the side band sample contains a 
contribution where the two mass assignments of the D° decay products are reversed. 
This contribution induces an asymmetry opposite to the real asymmetry. Its size 
(around 5% of the side band sample) is estimated from the Monte Carlo simulation 
and then subtracted from the measured background asymmetry, using the Standard 
Model values for the charm and b asymmetries. The correction shifts the measured D° 
background asymmetry by 0.004. The combinatorial background asymmetries for all 
channels, measured at the three centre of mass energies, are listed in Table |j. 

The charm asymmetry is extracted by means of an unbinned maximum likelihood 
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Table 4: Combinatorial background asymmetries, measured on data, at the three centre 
of mass energies, from the side band of the mass peaks. 



Decay channel 


peak — 2 GeV 


peak 


peak + 2 GeV 


D*+ 


— > 7T + K 7T + 


0.001± 0.045 


0.0050±0.0092 


-0.013± 0.036 


D*+ 


— > 7T + K~7T + 7T° 


0.016± 0.024 


0.0063±0.0050 


0.045± 0.019 


D° - 


-> K~7T + 


-0.105± 0.060 


-0.005± 0.013 


-0.015± 0.048 


D*+ 


— > 7r+K~7r + 7r + 7r~ 


0.021± 0.014 


0.0035±0.0029 


0.011±0.011 


D*+ 


7T s +K-7r+(7T ) 


0.036± 0.028 


-0.0026±0.0057 


0.006± 0.021 


D+ - 


-> K~7T + 7r + 


0.016± 0.024 


0.0015±0.0052 


-0.004± 0.020 



Table 5: Sources of systematic errors on the measured charm forward-backward 
asymmetry. The total errors are obtained summing in quadrature the relative 
contributions. 



Source 


AA C FB (%) 






peak - 2 GeV 


peak 




peak + 2 GeV 


Fraction of D mesons / s ; g 


± 0.15 


± 0.10 




± 0.22 


Charm fraction f c 


± 0.10 


± 0.02 




± 0.12 


b mixing 


± 0.09 


± 0.17 




± 0.18 


Comb. back, asymmetry 


± 0.98 


± 0.20 




± 0.75 


Gluon splitting 




± 0.03 




± 0.05 


TOTAL 


± 1.00 


± 0.28 


± 0.81 



fit giving the following results at the three different centre of mass energies: 

A c FB (y/s = 89.37 GeV) = (-1.0 ± 4.3 ± 1.0)% 
A c fb(Vs = 91.22 GeV) = (6.3 ± 0.9 ± 0.3)% 



A c fb(Vs = 92.96 GeV) = (11.0 ± 3.3 ± 0.8) 



/o 



The first error is statistical and the second arises from systematic uncertainties as 
listed in Table |5]. The angular distribution of the tagged Z — > cc events at the Z is 
shown in Figure |], after background subtraction and acceptance corrections. Figure [5] 
shows the measured asymmetries as a function of the centre of mass energy together 
with the predictions of the Standard Model. 

The pole asymmetry as defined in Eq. 2 is extracted from the measured asymmetries 
at the three energy points by expressing them as a single measurement at the Z 
mass and applying a correction for the effect of initial and final state radiation, QCD 
corrections and photon exchange and interference. As pointed out in Ref. ||, the 
theoretical estimate of the QCD correction has to be rescaled to take into account the 
bias from the experimental cuts. In particular, the requirement of high momentum 
j3(*) s rem oves events in which hard gluon emission occurred, substantially reducing the 
correction. As computed from Monte Carlo, the relative QCD correction after selection 
cuts are found to be consistent with zero within a 0.14% absolute uncertainty. The 
initial state radiation correction can also be biased by the selection cuts, which are less 
efficient at lower centre of mass energies. This correction, calculated to be 14.9% using 



the MIZA program [10], is lowered to (14.1 ± 0.1)% when the selection bias is taken 



6 



into account. The final relative correction to be applied to the measured asymmetry 
at the energy of the Z mass is (13.2 ± 0.2)%. 

Within the Standard Model, the measured asymmetry can be used to extract a value 
of sin 2 9$, taking into account the dependence of the b asymmetry on this quantity, 
yielding 

sin 2 Of = 0.2321 ± 0.0016 . 



5 Conclusion 

The forward-backward asymmetry in Z — > cc decays has been measured at three 
different energies, at the Z peak and off-peak at ± 2 GeV: 

A c FB (y/s = 89.37 GeV) = (-1.0 ±4.4)% 
Afb{Vs — 91-22 GeV) = (6.3 ±1.0)% 
A c fb(Vs = 92.96 GeV) = (11.0 ± 3.4)% . 

A Standard Model fit to the measured asymmetries yields sin 2 9$ = 0.2321±0.0016. 
These measurements are in agreement with the Standard Model predictions and with 



the other determinations of these parameters at LEP [11]. 
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Figure 1: Mass- difference distribution for candidates of the decay channel D* + — > 
7T+D° ; followed by D° -> K-tt+^; ; D° -> K-tt+ttY&J, D° -> K-tt+tt+tt- fcj and 
D° — > K _ 7r + (7r°)('d y ) . The full histogram are the data while the dots with error bars 
are the combinatorial background taken from Monte Carlo simulation. The error bars 
represent both the statistical error due to the limited Monte Carlo statistics and the 
systematic error in the normalisation. The dotted lines show the selected region. 
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Figure 2: Mass distribution for candidates of the decay channel D° — > K~7r + . T/ie 
dote wit/i error bars are data while the solid line represent the best fit to the distribution 
and the dashed line is the best fit to the background. 
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Figure 3: Mass distribution for candidates of the decay channel D + — > K~7r + 7r + . The 
dots with error bars are data while the solid line represent the best fit to the distribution 
and the dashed line is the best fit to the background. 
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Figure 4: Normalized angular distribution of the tagged Z — > cc events at y/s = 
91.2 GeV, corrected for acceptance and background subtracted. The dots with error 
bars are data while the dashed line is the result of the fit. 
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Fi gure 5: A c p B as a function of the centre of mass energy. The dots with error bars 
are the measured asymmetries. The curve is the prediction of the Standard Model with 
m t = 175 GeV/c 2 , m H = [90, 1000] GeV/c 2 and a s = 0.120. 
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